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Abstract. For determining distance and mass of pulsating stars a new, purely photometric method is described in 
which radial velocity observations are not needed. From multicolour photometry the variation of angular diameter 
is determined in a conventional way by using the surface brightness of the theoretical atmospheric models ATLAS 
of Kurucz (1997). As a function of phase the following two parameters are introduced in the Navier-Stokes 
equation: gravity (of the appropriate static ATLAS model), and angular diameter. Distance and mass are derived 
from phases of standstill. Conclusions are drawn on the hydrodynamic behaviour of the atmosphere. The new 
method is compared with the Baade-Wesselink (BW) method. As an example the RR Lyrae variable SU Dra is 
given: using UBVRc photometry, a distance of 647 ± 16pc, and a mass of 0.66 ± .03A4q were found. In addition 
to the radius change (4.58 — 5.51)-Rq radius undulations have been found with amplitude ~ 0.27?q and period 
P/5, which are synchronized with the main period P — 0"?66042. The present method has confirmed the distance 
value from the BW method; thus, in the problem of different distance scales from BW and other methods the 
short extragalactic distance scale has been bolstered by an independent argument. The mass is some 30 percent 
larger than the value from the equation of pulsation of van Albada & Baker (1971). 



Key words, stars: variables: RR Lyr - 
individual: SU Dra 



stars: fundamental parameters: distance, mass - stars: atmospheres - stars: 



1. Introduction 

The availability of ATLAS stellar atmospheric models 
(Kurucz, 1997) made it possible to account for the 
monochromatic flux of non- variable stars with a formerly 
unprecedented accuracy. As a function of the parameters 
effective temperature Tg, surface gravity g, chemical com- 
position, etc, large grids of models are now available for 
different physical assumptions on convection and turbu- 
lence. Ab initio knowledge or reconstruction of photomet- 
ric filter functions made it possible to compute absolute 
fluxes in a variety of photometric bands and to fit them to 
stellar magnitudes and colours by an appropriate choice of 
zero points. The comparison of observed colours of stars 
with computed colours of model atmospheres made it 
possible to determine the physical parameters of a star. 
These new possibilities opened the way to explore the 
atmospheric parameters of variable stars which pulsate 
in such a manner that quasi-static approximation (QSA, 
Buonaura et al. 1985) is valid. 

The synthetic colours UBVRIJK of new and more so- 
phisticated versions of the ATLAS models with increased 
microturbulence were used by Liu & Janes (1989, 1990) 
to determine fundamental parameters of RR Lyrae stars 
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using the Baade-Wesselink (BW) technique. There are, 
however, many problems with the BW method: a theoret- 
ical overview of these was given by Gautschy (1987) while 
practical questions and the application for RR Lyrae stars 
were discussed in detail by e.g. Liu & Janes (1990), Jones 
et al. (1992). 

The present paper will approach the determination of 
distance and mass of pulsating stars, especially RR Lyrae 
stars, from a modified point of view compared to a conven- 
tional the BW analysis. The synthetic colours and fluxes in 
different photometric bands of ATLAS models (Castelli, 
1999) will be used to determine the variation of stellar 
angular diameter by a similar, slightly extended method 
as in the BW analysis. However, the angular diameters 
and values of the effective gravity g,, will be introduced 
in the Navier-Stokes (N-S) equation averaged over the 
continuum-forming layers: the distance and mass of the 
star will be determined from the momentum balance of 
the variable atmosphere . In this manner the most time 
consuming radial velocity observations are not necessary, 
their problematic conversion to radius changes becomes 
unnecessary, and the fundamental parameters distance 
and mass will be obtained simultaneously by a procedure 
which is different from the BW technique. 
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Section 2 describes the method generally and its adap- 
tation to UBVRc photometry which is a compromise be- 
tween the number of available observations and the ef- 
ficiency in determining physical parameters of a pulsat- 
ing stellar atmosphere. Sect. 3 summarizes the light and 
colour curves of SU Dra (=BD -h67°652) from UBV{RI)c 
observations and the results from the method as an exam- 
ple. In Sect. 4 the results are discussed and compared with 
those from BW analysis. The conclusions are given in Sect. 
5. 

Spherical symmetry of the pulsation (excitation of 
radial modes) will be assumed. Cgs units will be used 
throughout; departures from them will be indicated, r, 6, (j) 
will denote spherical coordinates with their origin in the 
stellar centre of mass. 

2. The method 

Four steps are involved. 

1. The obscrv(xi magnitudes arc converted to fluxes at 
zero air mass and corrected for interstellar extinction. 

2. The observed coloiir indices arc converted to physical 
quantities of the theoretical atmospheric models: i.e. 
to Te, (/e, and surface flux in physical units for the 
different photometric bands including bolometric flux. 
In this step the bolometric flux is obtained from Tg 
by the Stefan-Boltzmann law while in the photometric 
bands the zero point of the magnitude scale must be 
determined to obtain the flux in physical units. 

3. From the combination of observed and model fluxes 
the variable angular diameter is obtained which is in- 
troduced in the N-S equation together with the values 

9e- 

4. At a flxed phase (time interval) the N-S equation is 
an algebraic equation for the unknown quantities dis- 
tance and mass. It must be written out for two phases 
at least, from which the unknowns are obtained by el- 
ementary operations. 

2.1. Selection of the photometric system 

A photometric system with colour indices C, is appro- 
priate; for the present method if the physical quantities 
are single-valued functions Te{Ci,C2), lgg{Ci,C2),.-. and 
they change slowly, i.e. the flnite differences \ATe/ ACi\, 
\Ag/ACi\, I A surface brightness/ ACj I are small for one of 
the indices Ci at least. E.g. in the ranges 6000 < Te < 
8000K, 2 < Ig 5e < 4 of an RR Lyrae star for UBV pho- 
tometry 

ATe/A(B - F)|c/-B=const « -10*K/mag, 
Alg5(e/A(J7 - B)|s-y=const » -lOdex/mag (1) 

which result in errors w 100 K, 0.1 dex for a typical error 
0.01 mag in the colour index. For BVRc photometry 

Alg.ge/A(S - V^)|y-ijc=const ~ 70dex/mag, 

Alg ge/A{V - i?c) |B-y=const « -70dex/mag (2) 



which make this photometry unflt to determine Qe- The 
final result is that Te, Igge and the fiux in different pho- 
tometric bands are found with reasonable accuracy from 
a pair of the colour indices. The conversion — > Tq etc 
can be made conveniently using the Kurucz tables (1997) 
in which the interstellar reddening is included. 

2.2. The photometric half angular diameter 

We list some useful relations for convenience. At a phase 
ip in the photometric band x the luminosity of the star is 

L^{^) = 47ri?g(^)^^(^) (3) 

where the radius i?o belongs to zero optical depth in the 
reference frame with origin in the stellar centre of mass, 

J',= dXS^iX) dTAVciBA(rAj (4) 

is the physical flux (Kurucz 1979) on the stellar surface 
in the photometric band x where Sx{X) is the filter func- 
tion defining the colour system. The monochromatic flux 
on the stellar surface was approximated in Eq. (4) by a 
sum, where ci,...,ce and TXj^,...,TXg are 0.1615, 0.1346, 
0.2973, 0.1872, 0.1906, 0.0288, and 0.038, 0.154, 0.335, 
0.793, 1.476, 3.89 respectively (Traving et al. 1966), Bx 
is the Planck function which accounts well for the source 
function if the stellar continuum is considered. The mean- 
ing of Eq. (4) is that the monochromatic flux originates 
roughly from six layers between Rt^kis.sb ^ r < RrxxiO.oia 
and the main contribution comes from the neighbourhood 
of i?rA~o.335- These (wavelength dependent) radii are de- 
fined e.g. by 

/'^iT;^=0.335 

rx = 0.335 = - / nx{r)dr, (5) 

jRo 

Kx is the monochromatic absorption coefficient. If Kx{r) 
is a smooth and moderately changing function of A we 
can assume that Ri = Rtx=().q3h- R = Rtx=o.335, R2 = 
Rtx=3-S9 are independent of wavelength (overlining de- 
notes averaging over A). If x represents bolometric or an 
optical band the Rosseland optical depth can be substi- 
tuted for Tx to estimate R, R2, which are now au- 
tomatically independent of the wavelength. Rq is the top 
of the atmosphere, Ri, R2 can be regarded as the top 
and bottom of the continuum-forming atmospheric region. 
(i?i - R2)/Ra, R/Ro vary with (p, R2 < R < Ri. These 
variations reflect the change of Ig (7e ((/?)• 

At distance d the observed stellar flux Ix will be 

Ix{v) = 10-^-/2.5^^(^^/4^^2 (g) 

where Ax is the interstellar extinction in magnitudes. The 
half angular diameter of the star is obtained for zero op- 
tical depth by 

M^) = M'P)ld = [10^==/^%(^)/.?.((/.)]i/2. (7) 
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Xx must be derived from observed magnitudes and com- 
pared with the tabulated value Tx of the theoretical mod- 
els which are found by a comparison of observed colour 
indices with those of computed models. 

The conventional way to determine i^oliy^) is to intro- 
duce A„ Ix{^) = ^olO-['».('^)+BC.(.p)]/2.5^ and Jfbol(^) = 

oT^i^^)) in Eq. (7) where a is the Stcfan-Boltzmann con- 
stant and the zero points of the magnitude scales nix-, 
BCa; are merged in the constant S'o (e.g. Liu & Janes 
1990, X = V), i.e. two parameters of theoretical models 
(BCy, Te) belonging to a pair of observed colour indices 
are needed. This procedure was applied in the present pa- 
per as well; the actual form was 



^0 = 0.636 X 10-(>'+BCv-2A^)/5j.-2 



in radians. To check Eq. (8) a single parameter determina- 
tion of '&Q was used by the following formulae in Eq. (7). 
An apparent magnitude X is equivalent to a physical flux 



Zx = 10-^+(^^^«--^)/2-565'''*5'^^crgsscc-i cm-2, 



(9) 



where Xvega = 0.03, 0.039 mag, t^''*^''' = 1.82, 1.04 \iX = 
V, Rc, respectively. In the band X the physical flux of a 
model is 



(10) 



where Xk is the tabulated magnitude of the Kurucz (1997) 
tables for U, B, V, Rc, Ic- The derivation of Eqs. (9), (10) 
is given in Appendix A. 

2.3. The momentum balance of a radially pulsating 
atmosphere 

For a non-rotating star the r component of the N-S equa- 
tion is 



dvr{r,e,(l>,t) dvr , a M ^\ 



GM 



1 dp{r,e,^,t) 



p(r, 6, (j), t) dr 



■a2{r,e,.p,t) (11) 



(Landau & Lifshitz 1980) where is the r component 
of the velocity, G is the gravitational constant, A4 is the 
stellar mass, p is the density, ai = [vg{r,6,(j),t)dvr/d6 + 
sin~^Ov^{r, 9, (f), t)dvr/d(p-Vg{r, 9, 0, t)-v'^{r, 9, 0, t)]/r is 
the rest of the terms from {v, V)^^. The sum of the terms 
from the molecular viscosity is 02. Setting the sound ve- 
locity (about 10^ cm/sec) for the components vg, as 
an upper limit, the order of magnitude of ai will be 10 
cm/sec^; therefore, except for phases of stationary flow 
(i.e. dvr/dt « 0) this term is negligible in comparison with 
the explicitly written terms. Estimations are not avail- 
able for a2, however, it can be assumed to be negligible in 
comparison with the other terms, especially at phases of 
standstill (i.e. Vr ~ 0). 

For a static stratification the acceleration 



l&p 
p dr 



(12) 



is constant within few percents in the domain R2 < 
r < Ri. On the right hand side of Eq. (11) the second 
term can be approximated by a constant ge which is ob- 
tained from the colour-colour diagram. To take all dy- 
namical corrections into account the term ad{r,9,4',t) = 

ai + Vrdvr / dr -\ will be introduced; '• • •' represents the 

terms which were not listed here, e.g. a dynamical correc- 
tion in Eq. (12). However, the spatial average gj, of aa is 
expected to be small in the phases of standstill in com- 
parison to the acceleration dVr/dt. 

Perfect spherical symmetry of the velocity field will be 
assumed; the terms describing eventual non-radial modes 
are assumed to be zero: 



(8) vr{r,t)=vi'^\t) + ^ 



r=R 



(^R-r) + ... = R{t) + ---,{lZ) 



the homogeneous velocity i? is a good first approximation, 
the next term would account for the compression or ex- 
pansion of the atmosphere, i.e. for a change of R\ — R2. 

By introducing Eqs. (12), (13) in Eq. (11) and averaging 
over the interval R2 ^ r < Ri (i.e. multiplying by r^sin^, 
dividing by 47r(-Rf — i?2)/3, integrating over r,9,<p) the 
equation 

dVr{R,t) GM 



+ ge{t) + 9d{R,t) 



(14) 



dt 

is obtained. 

Now all terms are neglected in Eq. (13) except for R 
and an apparent angular change 'da = Ra/d is introduced 
such that 

i?0(t) = R{t) + Ra{t) = m) + ^a(t)]d, (15) 

Ra <^ R accounts for the angular change by opacity as a 

consequence of the variation of g^. In a turning point of 
•&0, i-e. at standstill of the atmosphere, i?o ~ 0, ~ 0, 
\dvr/dt\ ^ \vrdvr/dr\, thus, Eq. (14) takes the form 

M 

[Mt) - Mt)]d = + 5e(i) + gd{t) (16) 

and gd is expected to be small. 

The functions '&o{t) and ge{t) can be determined from 
the obs{^rv(xl colour indices and brightness, ?9o can be ob- 
tained by numerical differentiation. Ra can be estimated 
by interpolation from the Kurucz tables (1979); compared 
to R{) it is negligible if \gge > 2.5. Ra can be estimated by 
numerical differentiation; it is negligible during the pulsa- 
tion except for the phases in which g^ changes strongly. 
From continuity of do{t) it follows that there exist two 
turning points at least: t\, t2- Writing Eq. (16) for t = t\ 
and t2, neglecting gd, and solving the two equations leads 
to: 



d = 



ge{tl) + Raih) - [ge{t2) + Ra{t2)mt2)l'd{h)f 

Mh) - Mh)['^{t2)/'»{t^)Y 



M = [ge{t) - ^{t)d]i)^(t)d'^/G. 



,(17) 

(18) 



Since i?a is a small correction to ^0 Eqs. (17) and (18) 
were solved while neglecting the apparent change of radius 
{Ra = 0, Ra = 0, i.e. i9 = '!?o) and taking into account 

Ra ^0, Ra^ 0, l^a = Ra/d. 
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Table 1. UBV{RI)c magnitudes of SU Dra on HJD = 
2452054, segment k = 52. The fraction of HJD is given for 

V , the epochs of B, U, Rc, Ic are obtained by adding —.0007, 
-.0018, .0005, .0010 to the epoch of V. 
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3. Distance and mass of SU Dra, some physical 
parameters of its atmosphere 

As an example some parameters of the RR Lyrae vari- 
able SU Dra will be determined. This star has a stable 
light curve without Blazhko effect; however, the analy- 
sis of the photometric observations loft some ambiguity 
whether secondary variations of some 10 mmag do ex- 
ist (Barcza 2002, Paper I). SU Dra is a metal-deficient 
Population II star: [M] = — 1.6dex. Its galactic coordi- 
nates are / = 133.45, h = -1-48.27 i.e. it lies high above 
the disk, the interstellar reddening is small in that direc- 
tion: E{B -V) = 0.015 mag (Liu & Janes 1990). Its pho- 
tometric behaviour was described in Paper I; from this 
observational material the segments k = 25 — 51 were 
used to compile an average UBV{RI)c hght curve because 
these observations were performed using two comparison 
stars. To improve the phase coverage a new segment k = 
52 was added containing 19 UBV{RI)c observations on 
JD = 2452054 with the Wright camera attached to the Im 
Ritchey-Chretien telescope of the Konkoly Observatory, 
Table 1 reports these results. Technical details of the ob- 
servations and reduction were the same as in segments 
A; = 42 - 51 (Paper I). 

3.1. Observed UBV{RI)c light and colour curves 

The observed magnitudes U, B, Rq were linearly interpo- 
lated to the epoch of the V observations to avoid system- 
atic errors of colour indices U — B, B — V, V — Rc in 
the rising branch. A comparison with the colour indices 
of Liu & Janes (1989) indicated small systematic differ- 





- ^ 
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Fig. 1. Light curve V and colour curves of SU Dra. Dots: seg- 
ments fe = 25 — 29, 42 — 51 from Paper I and k = 52, all are 
shifted by the amount indicated in the text. Crosses: segments 
fe = 30 — 41 from Liu & Janes (1989). Lines: fitted curves to 
segments k = 25 — 52. 



ences between the colour systems at Konkoly Observatory 
and Kitt Peak. Since it is of primary importance to have 
the same filter functions in the observations and com- 
putation of synthetic colours of the atmospheric mod- 
els, the colours of the Konkoly observations have been 
shifted by A{B - V) = -1-0.029, A{U - B) = +0.026, 
A{V — Rc) = +0.019 mag. By these shifts the brightness 
V and the colour indices of one of the comparison stars, 
BD -|-67°708, became identical in the Konkoly and Kitt 
Peak colour systems as is obvious from Table 3 of Paper 
I, and the systematic differences of the Konkoly and Kitt 
Peak curves disappeared. Fig. 1 is a plot of the observed 
magnitudes V and the colour indices, the phase was com- 
puted according to Paper I taking into account the secular 
change of period and phase noise. The curves drawn in 
Fig. 1 were obtained by high-order spline functions, they 
were fitted to the shifted colour cmves. These curves are 
the basis for deriving Te{(p), Iggd^), J-x{^), ^xif) with 
X = V, Rc^Cv , and '&o{(p) according to Sect. 2. 

The colour curves V — Ic were determined as well; 
however, they were omitted from the analysis because the 
CCD realization of Ic was not free of systematic errors; 
these originated from the red cut off which is considerably 
different from that of a GaAs cathode and of the improved 
5/c(A) of Bcssell (1990). 

The fitted curves show definite structure around (p ~ 
0.4, 0.7 in addition to the ascending branch which is their 
main feature. If only 60-70 points are observed (e.g. Liu 
& Janes 1989) the fine structure of the light and colour 
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Fig. 2. Panel (a): lines: from bottom to top the isograv- 
ity curves Igg = 2, 3, 4, dashed lines: from left to right the 
isotherm curves Te = 8000, 7000, 6000 K. The curves are in- 
terpolated to [M] = -1.6, E{B -V) = 0.015. Dots: U - 
B,B — V loop of SU Dra. Characteristic points of the loop 
ip = 0.015 and 1.015,0.3,0.42,0.56,0.81,0.897,0.92,0.947 are 
marked by different symbols which are plotted in Panels (b)- 
(d) as well. Panels (b), (d): the interpolated Igge, values. 
Panel (c): -do from Eq. (7), Eq. (9), Eq. (10) with V, Rc and 
from Eq. (8). The symbols of the loop indicate '&o{ip) from Eq. 
(8). 



curves remains hidden. It could be revealed by the 674 
observed points in the present study. 

3.2. Effective temperature, gravity, angular diameter 

The synthetic colours from the Kurucz (1997) tables 
(Castelli 1999) were interpolated linearly to [M] = — 1.6 
dcx and E{B — V) = 0.015. Te{(p), Iggeif), tabular mag- 
nitudes V, Rc, BCv were interpolated to a pair of colour 
indices {U-B,B-V), {U-B,V-Rc), {U-V,V-Rc). 
V, Rc were converted to surface flux J^y, ^Rc by Eq. 
(10). 

Panel (a) of Fig. 2 shows the theoretical isogravity 

curves \gg{U — B , B — V) = 2, 3, 4 for different values of 
Te, the isotherm curves Te{U - B, B -V) = Q,7,8 x lO^K 
for diff'erent values of Ig^, and the observed colour loop 
{U — B, B --V) as a function of (p. Panels (b)-(d) are a plot 
of IgQeif), 'i^oif), Te{(f) for the colour loop, t^olv) was de- 
termined from Eq. (8) and from {Iv,^v), (Irc^^Rc) ^y 
Eqs. (7), (9), (10): it is remarkable that the curves arc con- 
gruent within 2 percent. This finding demonstrates that 
in the range Ig^e, Te of an RR Lyrae star the simultaneous 
use of B — V and U—B estimates well. Therefore, the in- 
frared photometry can be dropped; it is more complicated 



from an observational point of view and was expected to 
give Te correctly by observing the slope of the continuum 
over a larger wavelength interval (Liu & Janes 1990). 

The procedure was performed with the colour-colour 
curves [U - B,V - Rc), (U - V,V - Rc). The con- 
gruence of the curves from Eqs. (8) and (7), (9), (10) 
was within 2 percent confirming that can be deter- 
mined well by simultaneous use of U — B oi U — V and 
V — Rc. The curves \gge{ip) etc. were of similar form with 
those from {U — B,B — V), the differences did not ex- 
ceed O.ldex, 100 K in the descending branch while in 
the ascending branch upper limits of 0.4dcx, 200 K were 
found for the differences. The curves '&o{'fi) were coinci- 
dent to within 10 percent. The differences were attributed 
to non-perfect filter functions, and to the fitting proce- 
dure when the averaged colour curves were constructed. 
Among Mv) from {U-B,B-V) and {U-B,V- Rc), 
{U — V,V — Rc) the largest differences were in the phase 
interval 0.87 < (p < 1; however, the cause for this seems 
to be that the curves based on B — V and V — Rc respec- 
tively are asynchronous. The differences could have been 
removed by a phase shift; nevertheless, averaging seemed 
more appropriate. Therefore, at the points Lp the values 
Iggei'p), Te{(p), ^oif) were averaged from the three com- 
binations of colour-colour indices of UBVRc- The result 
is plotted in Fig. 3. 

The values of lgge{(p), Te{(p) agree well with those 
Siegel (1982) and Liu & Janes (1990) determined from 
Stromgren and UBVRIJK photometry respectively. The 
averaged values arc \gge{ip) — 2.41 and Te{'p) = 6490 K 
(this study) while 2.69, 6400 K, and 2.72, 6433 K are given 
by Siegel (1982) and Liu & Janes (1990) respectively. 

The agreement between the present t?o(¥') curves and 
those of Liu & Janes (1990) is satisfactory. The main 
difference lies in the undulation of the present curve 
amounting to some 5 percent which was provided by ten 
times more photometric observations. The wavy struc- 
ture of t^oCv) is remarkable: at about w 0.45, 0.7 local 
minima of ^?o(<^) are clearly visible. This suggests that 
a sub-oscillation exists which is synchronized with the 
main pulsation. In the phase intervals 0.8 < (p < 0.9, 
0.1 < p < 0.4 the change of -da and eventual dynamical 
effects from the strong shock hamper us in getting a clear 
picture of the sub-oscillation; however, a periodic undula- 
tion with « P/5 is clearly visible in i9o{p) and V{p) where 
P = 5.71 X 10'' sec is the period of SU Dra. To check the 
reality of the undulations, the 674 V observations were 
analysed by the string length minimization described in 
Paper I since they can be regarded as the most reliable 
because of the use of two comparison stars. The period 
P/5 has been found only with a very shallow minimum of 
the string length. 

The averaged curve doif) was differentiated by mid- 
point formulae with a phase span of Ap = 0.02(?a 20 
minutes) which was found to be an appropriate compro- 
mise between the noise and the accuracy of i?o- This choice 
is in agreement with the present simplified hydrodynamic 
model: a description of more rapid structures cannot be 
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Fig. 3. Dots: averaged curves '3q{ip), Te(ip), lggc{f) from the 
loops {U-B,B-V), {U-B,V-Rc), (U-V,V-Rc). &o was 
derived from the averaged ^oif) by numerical differentiation. 
Crosses: variation of i?o('/') given by Liu & Janes (1990) from 
infrared photometry using a formula of type Eq. (8). 

expected. The angular acceleration ifo is plotted in the 
lower left panel of Fig. 3. 

To estimate the apparent radius change Rai'p) the 
geometrical depth of tross = 0.335 was interpolated to 
Iggei^p), Te{(p) from the Kurucz (1979) tables. Roif) = 
1.99 X lO^^'&oi^) km, R{ip), and Ri{tp) - R2{(p) as error 
bars are plotted in Fig. 4a. From the point of view of ap- 
parent diameter change and acceleration the critical phase 
intervals are 0.8 < ip < 0.9, 0.1 < (p < 0.4. The break of 
the slope of "ffoif) at (/? « 0.25 can be interpreted by the 
peak of Ra which is superimposed on the true atmospheric 
motion Ro{<f)- Ra was obtained from differentiating poly- 
nomial fits to Ra. Three points are noteworthy: 

* The secondary peak of '!?o at « 0.55 coincides just 

with the fall of the infrared brightness K which begins 
at phase delay ~ 0.6 from the visual maximum (Liu & 
Janes 1989). 

* Double standstill exists at minimal radius with differ- 
ences Aifi = 0.031, A-do = t?o(0.947) - i9o(0.978) ?a 
10-12 (i g 2 X 10*km). 

* The most prominent features of do are produced by a 
strong variation of Ra = Ro — R- 



3.3. Distance and mass 

At both turning phases the fine structure of t?o(</') shows 
two standstills with roughly equal i^o, the averaged val- 
ues of the quantities for Eqs. (17), (18) arc listed in Table 
2. Selection criteria were for these phases that the atmo- 



Fig. 4. Physical quantities characterizing the pulsation of SU 
Dra, d = 1.99 x 10^'' km was used to convert the angular quan- 
tities. Panel (a): dots: Ro{ip), crosses: R{ip) which belong to 
T = and TRoss = 0.335 respectively. The error bars indicate 
Ri{ip) — R-zi^p). Panel (b): dots: 'ff^d, crosses: from CORAVEL 
radial velocities (Liu & Janes 1989), triangles and circles: from 
radial velocities of metallic lines and H-,. (Oke et al. 1962) 
1.32(— ?;radiai(¥') — 166.9) km/sec is plotted. Panel (c): dots: 
V - 9.2 crosses: Mboi = -2.5lglL{(p)/ Lq] 4-4.75. Panel (d): 
dots: ??od, line: -GM/R^ip), crosses: Oe- 

sphere must be in standstill to have the smallest gd in Eq. 
(14), the effective gravity must be constant to avoid strong 
variation of Ra. and the interval of the averaging must 
be long enough to smooth out scatter and eventual rapid 
variations for which the present hydrodynamic model is 
oversimplified. The errors of Ig^e and t9o arc less than 3 
percent; the errors of 'do and Ra dominate the error of 
d, M. In the phase intervals A and B the atmosphere is 
falling freely, in C and D it is in the state of maximal 
compression by the shock(s?) from the deeper layers. 

In the intervals A,B,C,D it can be assumed that the 
terms of gd with a factor v are negligible, i.e. ge of the 
static atmosphere accounts well for the outward accelera- 
tion via Eq. (14). The values ge etc of the pairs AC, AD, 
EC, BD, CD were introduced in Eqs. (17) and (18), d 
and A4 are listed in Table 3, columns 2,3 and 4,5 refer to 
Ra = 0, Ra = and Ra, Ra from Table 2, respectively. 
The pair AB was omitted because the errors in ge, ^o, 
Ra made d and A4 unrealistic. After averaging columns 
2,3 d = (564 ± 31)pc, M = (0.51 ± .06)^© were found 
as a first approximation; the errors represent 68 percent 
confidence level. 

Columns 4,5 are more surprising: they suggest 

d= (647±16)pc 

M = (0.66 ± .03)7^0. 
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Table 2. The averaged effective gravity, apparent acceleration 
of the continuum- forming layers, angular diameter, and accel- 
eration for the phase interval of standstills. The units are 
m/sec^, m/sec^, lO^'km, 10~^°radian, 10~^® radian/sec^. 









9e 


Ra 


Ra 


1?0 


i9o 


A: 


.305- 


.315 


1.12 


2.27 : 


16.9 


1.93 


-4.75 ± .56 


B: 


.545- 


555 


1.58 


-0.25 


11.3 


1.88 


-2.67 ± .45 


C: 


.941- 


954 


50.1 


3.45 


0.49 


1.61 


22.6 ± 2.5 


D: 


.973- 


.983 


31.6 


4.02 


0.99 


1.60 


13.6 ± 1.3 



: uncertain value 



Table 3. Distances and masses from Eqs. (17), (18) without 
and with correction for Ra, in units of pc and Mq respectively. 





d 


I 

M 


d 


11 

M 


AC 


535 


0.68 


555 


0.85 


AD 


476 


0.49 


512 


0.68 


EC 


593 


0.57 


647 


0.66 


BD 


554 


0.48 


647 


0.66 


CD 


665 


0.31 


647 


0.66 



I: Ra = Q, Ra = 

II: Ra and Ra from Table 2, ■da = Ra/MYpc 



The values using A were excluded from these averaged d, 
M on the following arguments. The atmosphere is very 
extended in the interval 0.15 < cp < 0.4 i.e. its optical 
scale height increases to i?o — -^2 ~ 3 x lO^km « 0.1i?o at 
(f = 0.225, plane parallel approximation (used in ATLAS 
models for radiative transfer) is questionable (Baschek et 
al. 1992). The apparent acceleration exceeds the effective 
gravity {Ra ~ ^9e), ~ is presumably not satisfied, 
and it begins just a strong apparent (linear rather than 
accelerating?) contraction; furthermore, Ra can be esti- 
mated with large uncertainty. In Table 2 Ra{v = 0.31) = 
2.27m/sec^ is probably an overestimation, e.g. an ad hoc 
correction to Ra{f — 0.31) = O.lOm/sec^ would result 
in d = 587, 559pc and M = O.81,O.7OAl0 which arc 
somewhat closer to the values from combining B,C,D. 
Therefore, the data from AC, AD have a very low weight; 
however, to give a formal standard error to the unequivo- 
cal d and A4 from BC,BD,CD they were used but omitted 
in the averaging. 

3.4. Luminosity, radius, velocity, acceleration 

The distance modulus is 

TO - M = (9.05 ± .05) mag. 

Using Roi'p) and T(.{(p) the intensity averaged luminosity 
is (L) = (1.68 ± .08) X lO^^ergs/sec = (42.4±2.1)Lq; the 
errors were derived from the error of d. The mean radius is 
Rq = 3.54x lO^km. The intensity and magnitude averaged 



visual absolute magnitudes are (M™*) = 0.74 ± .05 and 
{My''^) = 0.78 ± .05 respectively. 

Fig. 4b-d arc plots of different physical quantities de- 
rived with d= 1.99 X 10^^ km, the data were averaged for 
intervals Af = 0.01. Note that the velocity '&o{'P)d shows 
a correlation with p[v-y — Wradiai (</')] only if it is averaged 
over a longer time, {v^ = — 166.9km/sec (Liu & Janes 
1990) is the barycentric velocity, p = 1.32 is the conver- 
sion factor.) The phase difference between bolometric and 
visual maximum isA(p«0.015. 

4. Discussion 

4.1. Appraisal of the method, comparison with the 
BW technique 

The surface brightness method to determine stellar an- 
gular diameters is known in two variants. One of them 
(Barnes et al. 1978) uses the empirically determined sur- 
face brightness of non- variable stars and applies it to vari- 
ables. The other variant is boimd to static theoretical at- 
mospheric models which are the sources of surface bright- 
ness and effective gravity. When applying them to variable 
stars a reliable surface brightness can be determined in the 
phases where a strong deviation from QSA is not expected. 
In the present paper the theoretical variant was used: the 
half angular diameter and the effective gravity were ob- 
tained from ATLAS atmospheric models as a function of 
phase. Concerning RR Lyrae stars arguments for the ap- 
plicability of QSA were discussed in detail by Buonaura 
et al. (1985). Step 1 is identical in the present and BW 
methods. Step 2 is nearly the same in both methods, the 
difference is only that here the values -Sq from effective 
temperature and bolometric correction to V are checked 
by a direct comparison of fluxes in V and Rq- 

Steps 3 and 4 arc principally different from the BW 
technique: all the data '&o{'p), deif) are used in a dynam- 
ical equation, i.e. in the r component of the N-S equation 
averaged over the continuum- forming layers. Some terms 
could be neglected and the pressure-density stratification 
of atmospheric models with QSA was introduced. An es- 
sential assumption is that the fimction = .ge(Ci,...) 
gives the uniform outward acceleration in i?i < r < R2, 
and the sum g^. plus gravity determines the motion of the 
atmosphere at R. A term g^ was defined to account for 
the dynamical corrections to the stratification of QSA. In 
the range Tg, Igg^. of RR Lyrae stars the ATLAS models 
are convective. Its presumably small effect on the momen- 
tum balance of the pulsation could be incorporated in g^; 
however, it was neglected by the use of Eq. (12). In stand- 
stills of the atmosphere gd is expected to be negligible and 
the averaged N-S equation is simplified to algebraic equa- 
tions for the unknown quantities distance and mass. In 
the present method the mass plays an inherent role; it is 
determined simultaneously with the distance. 

In the BW method the mass of a pulsating star plays 
an indirect role: in an equation of type Eq. (14) with 
ga = an assumed mass is introduced to derive ge{<f) 
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by differentiating the curve — pi'radiai(¥') numerically. This 
effective gravity is then used to select the atmospheric 
models from which the effective temperature is determined 
by a relation of type Te(Ci), and the variation of angular 
diameter is determined. The kinematic equation 

= / P^T - '?^radial(v)]f'dV (19) 

is solved for d by various techniques (e.g. Liu & Janes 1990, 
Cohen 1992, Jones et al. 1992), in the interval [(^i, (/?2] the 
atmospheres must be free of shocks. The primary error 
sources are and p; it is obvious that their uncertainty 
can result in a considerable error in d (Sabbey et al. 1995, 
Fernley 1994). The undulations reported in the present pa- 
per have a negligible effect on d if |i^2 — '/'il is sufficiently 
large. Another small correction comes from the fact that 
R{f) would be more appropriate than Ro{ip) if the ra- 
dius change is determined by using radial velocities from 
CORAVEL (e.g. Liu & Janes 1989) or from weak metallic 
lines (e.g. Oke et al. 1962). 

In the present method i'radiai(<p), P, and an eventual 
arbitrary phase shift of the radius change from photome- 
try and radial velocity (e.g. Aip = 0.02 for SU Dra, Siegel 
1982) do not enter at all, the observation of radial veloc- 
ities is not needed. However, a colour must observed, e.g. 
U in Johnson system, which is a good indicator for grav- 
ity. In the BW method gd ^ can result in a negligible 
error, in the present method at selected phases gd — 0, 
i.e. the unconditional use of Eq. (12), and the simplifi- 
cation of the velocity field to Eq. (13) with a single con- 
stant term can be assumed as a source of more uncertainty 
than in the BW method which involves the assumption 
Vr{r,t) = R{t) as well. Problems arising from the differ- 
ence S^{X) - S°{X) ^ were discussed by Bessell (1990) 
where S°{X) represents the filter functions of an observa- 
tory. They are common both in the present method and in 
the BW method. This eventual systematic error is phase 
dependent. 

4.2. The results for SU Dra 

The distance derived in the present study agrees vc;ry well 
with its value d = 640 pc from the BW analysis (Liu 
& Janes 1990) if the subtle changes of Raif) are taken 
into account. It is in accordance with the low absolute 
brightness of an RR Lyrae star if it is determined from 
a direct analysis. The higher absolute brightness of an 
RR Lyrae star from cluster properties is not bolstered 
by the present value of d. Concerning the short and long 
extragalactic distance scale (Gratton 1998) the discrep- 
ancy has been strengthened and an additional indepen- 
dent argument has been given for the short scale. If the 
change of Ra{^) is neglected the discrepancy gets worse 
hy d = 564 ± 31pc. Furthermore, Eq. (19) has been ver- 
ified: p = 1.32 and Vj = —166.9 km/sec (Liu & Janes 



1990) must essentially be correct since e.g. for (fi2 = 0.8, 
(fii = 0.4 an error Av-y = 1 km/sec would result in an error 

Ad/d = [(^2 - ^i)pPAvj]/[R{^2) - R{^i)] ~ 0.04 (20) 

i.e. roughly —0.09 mag in My- (It must be noted that 
Oke et al. (1962) argue for v-y = —161 km/sec). In the 
BW method the secondary error sources undulation and 
use of R for Rq result in \Ad/d\ < 0.005. 

Liu & Janes (1990) assumed M = 0.55 — 0.65 to get 
Ig^e for the interval 0.4 < if < 0.8; they remarked that 
their d was insensitive to the actual value (0.6A^q). The 
basic pulsation equation of van Albada & Baker (1971) 
gives A^suDra = O.47A^0, from which the present value 
0.66 ± .03A^Q differs significantly, however, its accep- 
tance must be supported by a smaller uncertainty in input 
physics. 

The observed — pwradiai (</') (Liu & Janes 1989, Oke et 
al. 1962) does not show the undulation of dod, and the 
assumption of an eventually phase dependent conversion 
factor (Sabbey et al 1995) cannot bring the curves to 
congruence. For a few variables in M5 and M92 a simi- 
lar undulation of i?o(v) was shown graphically by Cohen 
(1992); however, it was not discussed. On a speculation 
level the most probable cause of this inconsistency is 
the non-negligible velocity gradient coupled with undu- 
lation of do- The radial velocities (especially those from 
CORAVEL) refiect the averaged velocity from the deep 
layers tross ~ 0.335 while "dod is the averaged velocity 
over the layers « t < tross = 0.038. At phase cp w 0.75 
a weak hump is discernible on the CORAVEL curve which 
is correlated with the local maximum t?o(0.78). It is, how- 
ever, puzzling why there is such a loose correlation of 

^oi(p)d with p[Vj - Uradial(<<3)]. 

5. Conclusions 

In the frame of quasi-static approximation for the atmo- 
sphere a new method has been described to determine dis- 
tance and mass of pulsating stars. In elaborating the pho- 
tometry there have been minor improvements compared 
to the conventional way in the BW analysis: the use of 
bolometric magnitudes from T^, minus the bolometric cor- 
rection was supplemented with a comparison of physical 
fluxes from models and observations in a single filter band. 
It was postulated that at the phases of standstill the static 
ATLAS models give the effective outward acceleration in 
the atmosphere correctly. The effective gravities, the half 
angular diameters and their second derivatives with re- 
spect to time were introduced in the r component of the 
Navier-Stokes equation from which we found distance and 
mass while neglecting dynamical corrections in the stand- 
stills of the atmosphere and making a small correction 
for the apparent radius changes. The method presented is 
purely photometric; radial velocity observations and their 
sophisticated interpretation are not needed. 

As an example the light and colour curves of SU Dra 
have been interpreted. The agreement of the photometric 
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angular diameters of different origin was within the ex- 
pected error of a few percent. The distance we found is 
in excellent agreement with that from the BW analysis. A 
higher mass was found than that from pulsation equations; 
however, it is well reconcilable with the general mass val- 
ues of RR Lyrae stars. An undulation of angular diameter 
has been found which can originate from a surface wave, 
restricted to the outermost parts of the atmosphere. It is 
synchronized with the main pulsation. 

The results are encouraging for an application to dis- 
tant RR Lyrae stars, e.g. for members of globular clus- 
ters or nearby galaxies. The necessary photometric ob- 
servations are simple and can be done for large samples 
of these faint stars by moderate efforts. Distance deriva- 
tions based on the BW method are available for some 30 
field and a few cluster RR Lyrae variables (Gratton et al. 
1997). These limited numbers could easily be increased by 
an order of magnitude if the present method were applied. 
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Table A.l. Magnitudes and physical fluxes of Vega in units 
10~®ergs sec~^ cm~^ 



X 


U 


B 


V 


Rc 


Ic 


magnitudes 


0.022 


0.027 


0.03 


0.039 


0.035 


-7-(Vega) 




32.9 


18.1 


10.5 


5.95 


,o(Vega) ^(Vega) 
''o •'' X 


21.6 


31.7 


18.2 


10.4 


5.87 



Appendix A: The zero point of the magnitudes 
UBV{RI)c of the Kurucz tables 

Tiig et al. (1977) measured the physical flux X^^'^^a) 
Vega at 90 wavelength in the interval 3295 A < A < 

9090 A (Flagstaff flux calibration for a Lyr); this is Ix 

for Eq. (7) in units ergssec"^ cm~^ A . The theoretical 
Eddington flux of Kurucz (1997, file fvega.pck belong- 
ing to Te = 9550K, Ig g = 3.95, [mctalHcity] = -0.5dcx, 
^^turb = 2km/s) was multiplied by 47r, converted to the 
same units and interpolated to the 90 wavelengths: this 
is for Eq. (7). Ax was set to = 0, ^^^"^'"'^ = 

(7.95 ± .02) X 10~^ was found which agrees well with the 
interferometric value (7.85 ± .17) x 10~^ (Hanbury Brown 
et al. 1974). Having BCy = -0.310, = 9550 K from 
the theoretical model and V = 0.03 mag from the obser- 
vations the constant 0.636 is obtained in Eq. (8). (The 
parameters differ slightly from the observational values of 
Code et al. (1976): BCy = -0.250, = 9660 K.) 

The observed flux X^^°sa) ^-^^ theoretical flux 

^(Vega)^|Vega) integrated with the filter functions for 
the system UBV{RI)c (BcsscU 1990) normalized to unity 
to obtain the physical fluxes belonging to the magnitudes 
of Vega. Table A.l shows the results. lu could not be com- 
puted from the Tiig ct al. (1977) calibration because of a 
lack of data in the ultraviolet. The difference of l''^'^^^^ 
and d'^^'^^^^ J^'^^^^'^^ originates from hues: J^^°^^^ was not 
measured at the strong lines and Balmer jump; it was in- 
terpolated linearly. 

For computing the observed flux of a star at zero air 
mass the lower row of Table A.l must be used, an apparent 
magnitude X is equivalent to a physical flux 

= 10(^Vega-X)/2.5^(Vega)^O'ega)g^gggg^_l ^^_2_ 

The monochromatic Eddington flux Hx of Vega 
(Kurucz 1997, fvega.pck) was integrated with the 
same normalized fllter functions UBV{RI)c of Bessell 
(1990); -2.51gHx = -23.835, -24.250, -23.645, -23.083, 
—22.370 were foinid for X = U, B,V, Rc, Ic respectively. 
On the other hand the tabular Vega surface magnitudes 
are Xk = -19.339, -19.755, -19.151, -18.590, -17.875 
if the improved filter functions of Bessell (1990) are used 
(Castelli 2002). The difference in these magnitudes is 
—4.495 i.e. the physical flux of the model is 
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where Xk is the tabulated magnitude U,B,V,Rc,Ic 
(Kurucz 1997). 



